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Abstract 



Hypercycles are information integration systems which are thought to overcome the information crisis of prebiotic evolution 
rj-J by ensuring the coexistence of several short templates. For imperfect template replication, we derive a simple expression for the 
p maximum number of distinct templates n m that can coexist in a hypercycle and show that it is a decreasing function of the length 
. "L of the templates. In the case of high replication accuracy we find that the product n m L tends to a constant value, limiting thus 
O the information content of the hypercycle. Template coexistence is achieved either as a stationary equilibrium (stable fixed point) 
or a stable periodic orbit in which the total concentration of functional templates is nonzero. For the hypercycle system studied here 
I we find numerical evidence that the existence of an unstable fixed point is a necessary condition for the presence of periodic orbits. 
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1. Introduction 

Most of the modern theoretical work on pr ebiotic evolu- 
tion was prompted by the seminal paper of lEigenl (|197lh 
which explored the fate of a population of competing 
macromolecules in an environment with limited resources. 
The main conclusion of Eigen's work is that the length 
of a replicating polymer (i.e., a RNA-like template) is 
limited by the replication accuracy per nucleotide, and 
so primordial replicators would have to replicate with 
implausible high accuracy in order to reach the length 
of today's RNA viruses (about 10 3 to 10 4 nucleotides). 
This finding, together with the observation that distinct 
templates cannot coex i st in t he competition-only scenario 



(jSwetina and SchusterL Il982l ). come to be known as the 



information crisis of prebiotic evolution. 

In fact, this information crisis would be resolved by the 
coexistence of several distinct, short templates, i.e. by the 
splitting of the information in short modules, similarly to 
the division of the genome in chromosomes found in many 
present-day organisms. In this case, the total information 
content of the template pool is the product of the num- 
ber of different templates and the maximum information 
coded per template (roughly the template length L), pro- 
vided the template types have the same concentration. But 
template coexistence can be achieved only if some sort 
of cooperation between templates is imposed a priori to 



the m olecular population. To this end Eigen and Schusteil 
19791) proposed a cyclic reaction scheme, termed hyper- 



cycle, in which each replicator would aid in the replica- 
tion of the next one, in a regulatory cycle closing on itself. 
An alternative proposal confines the templates in pack- 
ages or prebiotic vesicles which are deemed viable pro- 
vided it encloses a n umber n of distinct functional tem- 
plates dNiesert et ail 



Zintzaras et 



2112002; 



981; Szathmarv and Dcmeter, 1987 
Fontanari et al.l . l2006h . ~ 



Perhaps because hypercycles and vesicle models are more 
difficult to analyze than the naked-gene (quasispecies) sce- 
nario, the all-important question of whether these informa- 
tion integration systems exhibit a similar phenomenon as 
the error threshold of the quasispecies model was put off. 
Only recently a re-ex amination of a prot otypical package 



model - the model of Niesert et al. ( 198 lh - revealed that 
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package models in general suffer from the same malady as 
the quasispecies model: in the case of imperfect replica- 
tion an increase in the number n of distinct templates con- 
fined in the vesicle must be followed by a decrease of their 
lengths, otherwise the package becomes unviable. As a re- 
sult, the product nL (i.e., the total information content of 
the package) tends to a constant value that depends es- 
sentially on the spontaneous mutatio n rate per nucleotide 
( Silvestre and Fontanaril 120071 12008T ) . Our aim in this con- 
tribution is to investigate whether a similar restriction to 
the total amount of information in the pool of templates 
holds for hypercycle systems as well. 

The dynamics of hypercycle systems in the presence 



of a var iety of mutant types was extensively investi- 
gated by IStadler and Schuster (|l992fk lHappel and Stadler 
(Il998l). We refer the rea der t o bresch et al.l (Il980h: 



*A„ 



(3) 



Szathmarv and Demeterl (|1987t k ISardanves and Sole 
<|2007h for an emphasis on the destabilizing ef f ects o f 
mutant parasites and to iBoerliist and Hogewed (|199ll ): 
Cronhiort and Blomberg ( 1994 ) for an analysis of the ro- 
bustness conferred by spatial organ ization against those 
mutants. However, the formulation of Campos et al. (2000) 
in which the mutants form an error tail class seems more 
appropriate to study the error threshold phenomenon and 
so our analysis will build heavily on that paper. 

The sole motivation of this contribution is to show that 
the condition f or the viability of the hypercycle derived in 
Campos et al. (2000) is in fact valid for all n and not only 
for the regime where fixed points are stable, i.e., for n < 4. 
This is so because we found numerically that a necessary 
condition for the presence of stable periodic orbits is the ex- 
istence of an unstable fixed point. For the purpose of com- 
pleteness, in the following section we describe the model 
and re-derive the main results regarding the existence of 
fixed point (and hence of stable periodic solutions) for the 
hypercycle system. 



2. Model 

The hypercycle system we consider here is composed of 
n 'functional' elements I±, . . . ,I n and its error tail I e . The 
templates are capable of self-replication with productivity 
values Ai [i = 1, . . . , n) and A e . As usual, we introduce the 
kinetic constants Ki that measure the strength of the influ- 
ence of template Ii-\ on the growth promotion of template 
7j. The key ingredient in the modeling is that in both pro- 
cesses of growth of template li the probability of success 
is given by the parameter Q 6 [0, 1], so that an erroneous 
copy, which will then belong to the error tail, is produced 
with probability 1 — Q. Back- mutations from the error class 
to the functional class, as well as mutations between ele- 
ments of the functional class, are neglected. This formula- 
tion is equivalent to considering polynucleotides of length 
L — ► oo whose mutation probability per nucleotide u goes 
to such that the replication accuracy per genome is finite, 
i.e. exp {—Lu) — > Q. 

The concentrations Xi (i = 1, . . . , n.) of the hypercy- 
cle functional elements and the concentration x e of the 
error-tail evolve in tim e according to the kinetic equations 
(|Campos et alll2000l ) 



±i — Xi (AiQ + KiXi^iQ - $) i = l,...,n (1) 



and 



x e = x e (A e - $) + (1 - Q) ^ Xi (A, + KiXi^x) (2) 
where xq = x„ and 



is a dilution flux that keeps the total concentration con- 
stant, i.e., X^i=i x i + x c — 0. As usual, the dot denotes a 
time derivative. Henceforth we will assume that 



; = 1 



= 1. 



(4) 



i=l 



In accord with the usual assumption of package models that 
functional templates and parasites are selec t ively neutral 
(jNiesert et all 1 1 9 8 it ISilvestre and FontanariL l2008h we set 
Ai = A e = a, and Ki = K > for i = 1, . . . , n, resulting 
in the so-called symmetric hypercycle. Hence by measuring 
the kinetic constant K and the time t in units of a we can 
set a = 1 without loss of generality (except for a = 0, of 
course). 

3. Analysis of the steady state 

Here we focus only on the fixed-point solutions Xi > 
for i = 1, . . . , n. In this case, the condition ±2 = yields 
$ = Q + KQxi which, inserted in the equations £3 = . . . = 
x n = ii — 0, yields x± = x% = . . . = x n . Using these 
results in Eq. §3§ we find that X\ is given by the roots of 
the quadratic equation 

nKx\ - KQx x + 1 - Q = 0, (5) 

which has two real positive roots provided the condition 

KQ 2 - An (1 - Q) > (6) 

is satisfied. We note that K > is a necessary condition for 
the existence of hypercyclic solutions x\ = . . . = x n > 0, 
since in the absence of catalytic couplings among functional 
templates (i.e., K — 0) Eq. ([5]) has no solution. 

The analysis of the roots of Eq. ((5]) and the numerical 
evaluation of the Jacobian eigenvalues indicate that the 
smaller root is always unstable whereas the larger root is 
(locally) stable for n < 4. In addition, the disordered fixed 
point Xi = 0, Wi and x e = 1 is always stable (|Campos et all 
120001) . 

It is well known that for n > 4 the steady state of 
the h ypercycle is charact e rized by stable periodic solu- 
tions (lEigen and Schusterl. 19791 Hofbauer and Sigmund . 
Il988t iHofbauer et all Il99lh where the concentrations Xi 
vary wildly reaching values dangerously close to zero, 
which would certainly doo m a finite-population system 
(jNuno and TarazonaLll994l ). Nevertheless, here we take an 
optimistic stance and assume that the hypercycle is an ac- 
ceptable information integrator even in the regime where 
it exhibits periodic solutions: the n templates do coexist 
in this regime after all, albeit with (very) distinct concen- 
trations. The issue is then to determine the region in the 
space of the parameters Q and K where the steady-state 
solution is such that X)"=i x n — 1 ~ x e 1% nonzero. 

The search for steady-state solutions with nonvanishing 
concentration of functional templates can be carried out 




Fig. 1. Time evolution of the total concentration of func- 
tional templates for n = 12, K = 1000 and (top to bottom) 
Q = 0.5,0.4,0.3,0.2,0.19639,0.1. The periodic solutions disappear 
at Q ~ 0.19639, the value at which the condition (0 for the exis- 
tence of a nontrivial equilibrium is violated. We have no proof for 
this remarkable coincidence. For lower values of replication accuracy 
the dynamics converges to the trivial fixed point x e = 1, which is 
always stable. 

straightforwardly through the numerical integration of the 
differential equations ([1])- ([2|) . The procedure is illustrated 
in Fig. Q] for n — 12: for a fixed value of K we decrease 
Q until we reach a threshold value below which the only 
steady-state solution is x e = 1. 

Remarkably, through this numerical analysis we find that 
the region of viability of the hypercycle (i.e., the region 
characterized by a nonzero concentration of functional tem- 
plates, regardless of whether the stable solution is a fixed 
point or a periodic orbit) is determined by the condition of 
existence of real fixed points, Eq. ((6|) . We are not aware of 
a mathematical proof for this result. For n — 2, index the- 
ory can be used to show that inside the region enclosed by 
a periodic orbit there must exist at least one fixed point so 
that the absenc e of fixed points precludes the exi stence of 
periodic orbits ( Guckenheimer and Holmes] . 1983 . pp. 51). 
It seems, however, that this r esult cannot be exten ded to 
n > 2. The theorem proved bv lHofbauer et alJ ( 1991 ) is not 
helpful either: it asserts that for a hypercycle system with 
no error tail (i.e., Q = 1) if the (unique) fixed point becomes 
unstable then there is a stable periodic orbit, whereas our 
conjecture is that if there is a stable periodic orbit then 
there must be at least one unstable fixed point. 

In summary, we find that the hypercycle system is viable 
provided the number of functional templates n satisfies the 
condition n < n m where 

KQ 2 

nm = W^QY (?) 

Since n m is a monotonously increasing function of Q 6 [0,1] 
and Q, in turn, is a monotonously decreasing function of 
L (recall that Q = exp(— uL)) we find that n m decreases 
with increasing L. Hence, all other things being equal if the 
number of functional templates n is increased then their 
lengths L must decrease accordingly so as to guarantee that 
n < n rn is fulfilled. This remark shows that an information 
preservation principle similar to that derived for package 
models holds for the hypercycle as well, which suggests a 



reconsideration of the whole approach based on the coexis- 
tence of distinct templates to address the information crisis 
of prebiotic evolution. 

4. Conclusion 

In our analysis we have opted for the choice of parame- 
ters that most favored the stability of the hypercycle. For 
instance, introduction of other interactions such as the cat- 
alytic promotion of the growth of the templates in the er- 
ror tail by functional templates - an assumption implicit 
in the package models - can only reduce the value of n m . 
(We have verified that the catalytic coupling between func- 
tional templates and the error tail does not produce any 
qualitatively new result.) Similarly, the parameter setting 
that corresponds to the elementary hypercycle in which 
Ai = 0,Vi but A e = 1 results also in a reduction of n m . 
Asymmetric hypercycles in which the templates have dis- 
tinct productivities Ai leads to internal competit ion and 
again to the decrease of n m ([Campos et al. , 2000). Hence 
Eq. must be seen as an upper bound to the maximum 
number of functional templates in the hypercycle. 

The effect that a change SL causes on n m is given by 



1-0 



SL. 



(8) 



In the case the functional templates have a high replication 
accuracy (i.e., 1 — Q « uL w 0) we find that the product 
n m L tends to a cons tant value, similarly to the findings fo r 
the package models ( Silvestre and Fontanari . 20071 . l2008h . 

We stress that Eq. (JT]) is worthful because its validity 
is not restricted to the regime where the nontrivial fixed 
point is stable: it holds also in the regime where the only 
stable solutions are periodic orbits. We provide only nu- 
merical evidence to support this remarkable finding which 
is based on the conjecture that the existence of an unsta- 
ble fixed point is a necessary condition for the presence of 
stable periodic orbits in the system of ordinary differen- 
tial equations ([T])-((2]). It would be very interesting to find 
a proof for this conjecture. 
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